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Pole trajectories connecting light and heavy scalar mesons, both broad 
resonances and quasi-bound states, are computed employing a simple coupled- 
channel model. Instead of varying the coupling constant as in previous work, 
quark and meson masses are continuously changed, so as to have one scalar 
meson evolve smoothly into another with different flavor (s). In particular, it 
is shown, among several other cases, how the still controversial A"q(800) turns 
into the established XcO, via the disputed Z) s (2317). Moreover, a xj, (3946) is 
predicted, which may correspond to the recently observed y(3943) resonance. 
These results lend further support to our unified dynamical picture of all 
scalar mesons, as unitarized qq states with important two-meson components. 
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After more than four decades, understanding the scalar mesons continues to pose 
serious difficulties to theorists as well as experimentalists. Still today, no consensus 
exists about the lightest and oldest structures in the scalar-meson sector, namely 
the a (/ (600) [1]) [2] and the re (^ *( 800 ) [1]) [3,4]. But also the discovery of 
the surprisingly light charmed scalar .D* (2317) [5], though giving a new boost to 
meson spectroscopy in general, has not contributed to the understanding of scalar 
mesons, as can be seen from the many different approaches to the £)* (2317) in 
the literature (see Ref. [6] for a representative, albeit not totally exhaustive, list 
of references). Here, we shall focus on a formalism which successfully describes all 
mesonic resonances, including the scalar mesons. 

In Ref. [7] it was shown that the .D* (2317) meson can be straightforwardly 
explained as a normal cs state, but strongly coupled to the nearby DK channel, 
which is responsible for its low mass. The framework for this calculation was a simple 
coupled-channel model, which had been employed previsously [8] to fit the S'-wave 
Kit phase shifts, and predict the now listed [1] 1^(800), besides reproducing the 
established 1^(1430). Furthermore, another charmed scalar meson was predicted 
in Ref. [7], i.e., a broad Dq resonance above the Dn threshold, somewhere in the 
energy region 2.1-2.3 GeV, which may correspond to the Dq (2300-2400) [1,9]. Also 
higher-mass D* and Dq resonances were foreseen [7], which have not been observed 
so far. 

The purpose of this Letter is to show the interconnection of the scalar mesons 
7^(800), (2300-2400), £» s * (2317) with one another, and also with the established 
Xco(3415) [1]. Moreover, the same interconnection will be demonstrated for the 
higher-mass recurrences of these scalars, thereby finding a candidate for the very 
recently observed F(3943) charmonium state [10]. For that purpose, we shall employ 
the above-mentioned coupled-channel model, but now for fixed, physical coupling, 
while quark and threshold masses will be varied. Thus, a continuous and smooth 
transition can be achieved from one scalar meson to another. Crucial here will be 
a mass scaling [11, 12] of the two parameters modeling the off-diagonal potential 
that couples the confined and decay channels. This way, these two parameters, 
identical to the ones used in Refs. [7,8, 11], suffice to reasonably describe a vast 
range of distinct scalar mesons. On the other hand, the confinement and quark- 
mass parameters are taken at their usual published values. 

Starting point is a simple, intuitive coupled-channel model, describing a confined 
qq system, coupled to one meson-meson channel accounting for the possibility of real 
or virtual decay via the 3 Pq mechanism. If the transition potential is taken to be a 
spherical delta function, the lxl inverse K matrix can be solved in closed form, 
reading [8] 
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where je,ne are spherical Bessel and Neumann functions, respectively, A is the 3 Pq 
coupling, a is the delta-shell radius, E n g c are the energies of the bare confinement 
spectrum, B n £ c are the corresponding weight factors, p is the on-shell relative mo- 
mentum in the two-meson channel, given by the kinematically relativistic expression 

4s p 2 =\ s - (M 1 + M 2 ) 2 1 \s - (Mi - M 2 ) 2 1 , (2) 
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and n is the ensuing relativistic reduced mass 
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As the present paper deals with scalar mesons, we have i = and £ c — 1 in Eq. (JTJ . 
Moreover, since only ground states and first radial excitations are considered here, we 
shall approximate the infinite sum in Eq. ((!} by two confinement-spectrum states 
plus one rest term, also sticking to the numerical values used in Refs. [7,8, 11], 
namely B m = 1.0, B\\ = 0.2, and B 2 i = E 2 \ = oo, with B 2 \/ E 2 \ = 1. As for the two 
confinement levels, we parametrize them by a harmonic oscillator [7, 11], i.e., 

E nl = (2n + 2.5) uj + m gi + m q2 , (4) 

where u = 0.190 GeV, m n = 0.406 GeV {n = u,d), m s = 0.508 GeV, and m c = 1.562 
GeV, as in previous work [7,8,11,13,14]. Finally, we assume a mass scaling of the 
parameters a and A given by [11, 12] 

a ij = constant , y/JIij = constant , (5) 

where the labels ij refer to a particular combination of quark flavors, and /ijj = 
m qi m q J {mq-^rriq.) is the corresponding reduced quark mass. This procedure ensures 
flavor invariance of our equations. Using then the values A ns = 0.75 GeV -3 / 2 and 
a ns = 3.2 GeV -1 from the fit to the S'-wave phase shifts in Ref. [8], we have fixed 
all our parameters, 1 which allows to show the predictive power of our approach. 
For the required input mesons masses, we take the isospin-averaged values [1] M n = 

0. 1373 GeV, M K = 0.4957 GeV, and Md = 1.867 GeV. 

Now we can compute pole trajectories in the complex energy plane for scalar reso- 
nances and (virtual) bound states, by searching the values of s for which cot So{p{s)) = 

1. However, instead of freely varying A as in previous work, we shall keep X ns fixed at 
its physical value of 0.75 GeV -3 / 2 , while changing instead one of the quark masses, 
as well as one of the meson masses in the decay channel. This way we can make one 
scalar meson turn into another. For instance, by letting 



m qi = m n + a(m c - m n ) 
M l = M w -a (M D - M n 



< a < 1 , (6) 



we smoothly change the k (ns) meson, coupling to the ttK channel, into the _D S * (2317) 
(cs), coupling to DK. The poles themselves are numerically found and checked with 
two independent methods, i.e., the MINUIT package of CERN [15], and MATHE- 
MATICA [16]. 

In Fig. Q one sees in one glimpse the nine trajectories 

a: K*(70A) - D*(2114) - £> s * (2327) -> Xco(3472), 

b: K*(1522) -> D*(2673) £> s * (2840) Xco(4015), (7) 

c: K *(1788) -> £> *(2841) -> £> a * (2923) -> x; o (3946), 



1 Note that we use here somewhat shifted confinement levels as compared to Ref. [8], namely 
the ones following from Eq. Q. This gives rise to a slightly lighter and broader k meson, and a 
heavier X o *(1430). 
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Figure 1: Scalar-meson pole trajectories in the complex energy plane. Dots represent 
predicted resonances or bound states. See text and Eqs. ()7I8|) for further details. 

where the numbers between parentheses are the real parts (in MeVs) of the respective 
resonance/bound-state poles, the corresponding imaginary parts being 



Before discussing the actual trajectories, a few remarks are due concerning the pre- 
cise values found for the pole positions. Clearly, for such a simple model without any 
fitting freedom, moreover covering a vast energy range, a very accurate reproduction 
of the masses and widths of all experimentally observed mesons cannot, and should 
not even be expected. In particular, the inclusion of only the lowest, dominant 
decay channel for each state will certainly reflect itself in one way or another. For 
instance, the much too small width of our #q(1788), which should correspond to 
the observed [17] Kq(1820), is probably owing to the neglect of the important Krj' 
channel. Furthermore, the somewhat too large mass of our x c0 (3472), as compared 
to the established [1] Xco(3415), may very well be due to the omission of vector- 
vector decay channels, which are relevant for charmonium ground states [18]. Note, 
however, that the latter discrepancy of 57 MeV is quite insignificant when compared 
to the huge coupled-channel shifts in charmonium recently found in Refs. [19,20]. 
Notwithstanding, a clear identification can be made of our broad _K"q(704), Dq(2114), 
and #5(1522) states with the listed [1] # *(800), D* (2300-2400), and #*(1430) res- 
onances, respectively. Here, one should also notice that we give the real parts of 
the pole positions of our resonances, which usually do not coincide with the ex- 
perimental masses resulting from Breit-Wigner fits when the widths are large. As 
for the remaining observed mesons, our £> s * (2327) is very close to the Z)* (2317), 
while our x' c0 (3946), with a width of about 60 MeV, seems a good candidate for the 



a: # *(-251) - £>S(-118) £> a * (0) Xco(0), 

b: K*(-95) -> D*(-183) - D s * (-220) - X c0 (-382), 

c: K*(-12) -> D*(-A5) -> D* (-57) -> x' c0 (-29). 



(8) 
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brand-new [10] charmonium state y(3943). Finally, we predict the two medium- 
broad charmed mesons -Dg(2841) an d -D* (2923), so far undetected, as well as the 
very broad states _Dq(2673), Z}* (2840), and Xco(4015), which will be extremely hard 
to observe at all. In any case, the predictions for the latter higher-mass states may 
change significantly when additional decay channels are taken into account. 

Turning now to the trajectories themselves, it is remarkable to observe that 
physical states with radically disparate widths can be continuously connected to 
one another in flavor. This is one of the reasons why scalar-meson spectroscopy is 
so intricate. Moreover, as we shall see below, states on the same mass trajectory 
can have different origins when viewed as qq states distorted by meson loops, which 
point will become clearer when we study Fig. El Anyway, the first radial excitations 
of the ns, cn, cs, and cc systems are all on the same trajectory in Fig. ^ i.e., the 
one connecting the Kq (1788) and x' c o(3946). 

In Fig. the lowest states for the various flavor combinations are displayed 



> 

0) 



Re(E) (GeV) 



1.00 
I 



2.00 
l_ 



0.65 



-100- 



-200- 



-300- 



Z)*(2114) 

0.9 



0.8 



0.7 



.0.9 
.0.85 



£> s * (2327) 




4.00 
l_ 



3.60 



d 



1.0 0.8 0.4 

Xco(3472) 



Figure 2: Pole trajectories of lowest states as a function of A. Numbers indicate 
reductions relative to the maximum value. 

again, but now also showing how the corresponding poles move when the coupling 
A is reduced from its fixed value. We see that the _K"o(704) and the Dq(2114) appear 
to find their origin in the continuum, corresponding to infinitely negative imaginary 
parts of their pole positions, while the _D S * (2327) and Xco(3472) are connected to the 
confinement spectrum, with poles on the real axis. This is quite surprising for the 
nearby pair Z)o(2114)-.D* (2327). However, even the physical -D s * (2317) itself can 
be either interpreted as a "confinement" state [11,21], or a "continuum" state [7], 
depending on tiny changes in e.g. the parameter a. What this figure also shows 
is an extremely delicate balance of coupling effects. With a small decrease of A, 
the Dq (2300-2400) and especially the k meson would become even broader and 
thus almost impossible to observe experimentally, while the _D S * (2317) would be a 
resonance or a virtual state instead of a quasi-bound state. 



5 



Finally, in Fig. |3]a direct transition of the ^(704) into the L>^(2327) is displayed, 
by letting m n — > m c , M n — > Mr, as in Eq. (jUJ), and moreover in a different fashion. 
Namely, instead of giving the pole positions in the complex energy plane, we now 
plot the corresponding real and imaginary parts as a function of the varying quark 
mass, as well as the proportionally changing threshold value. It is striking to see 
how the .K"q(704) resonance quickly turns into a virtual bound state, while its real 
part remains almost constant. Here, we probably see the kinematical Adler zero [21] 
at work, which rapidly moves away as one of the decay masses increases from M^, 
thus allowing the pole to approach the real axis. Then, the pole moves along the 
real virtual state, until it touches the threshold at about 1.76 GeV, after 

which it becomes a bound state. Notice again the tiny margin, at least on this scale, 
by which the _D S * (2327) is bound. 




Figure 3: Real and inaginary parts of Kq(704) pole turning into the _D S * (2327), as 
a function of varying quark mass. Straight dashed line stands for decay threshold 
(real). 

To conclude, in the present paper we have shown how several light and heavy 
scalar mesons can be linked to one another, by continuously varying some of the 
involved flavor and decay masses. This way, the common dynamical nature of the 
studied — and probably all — scalar mesons, as ordinary qq states but strongly 
distorted due to coupled channels, is further substantiated. Thus, tetraquarks and 
other exotic configurations are not needed in this context. Moreover, we deduce that 
labeling scalar mesons as qq states as opposed to dynamical meson-meson resonances 
makes no sense, in view of the tiny parameter variations needed to turn one kind 
of pole into another. Rather, scalar mesons should be considered nonperturbatively 
dressed qq systems, with large meson-meson components, no matter if one uses a 
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coupled-channel quark model [22] or e.g. the quark-level linear sigma model [23]. 
As a consequence, the spectroscopy of scalar mesons is much more complex than 
for ordinary mesons, with the total number of potentially observable states being 
different from the number of confined, bare qq states. 

In the course of this analysis, we have also found a candidate for the new charmo- 
nium state Y(3943) [10]. It is true that such a resonance, if indeed a scalar, should 
dominantly decay to DD, a mode which has not been observed yet. However, the 
reported decay Y(3943) — ► ujJ/^ is OZI-forbidden, so that it cannot account for 
the measured sizable width of V = 87 (±22 ± 26) MeV. 

We thank D. V. Bugg for enlightening discussions about scalar mesons in general, 
and the new charmonium state Y(3943) in particular. This work was supported in 
part by the Fundagao para a Ciencia e a Tecnologia of the Ministerio da Ciencia, 
Tecnologia e Ensino Superior of Portugal, under contract POCTI/FP/FNU/50328/2003 
and grant SFRH/BPD/9480/2002. 



References 

[1] S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B 592, 1 
(2004). 

[2] J. E. Augustin et al. [DM2 Collaboration], Nucl. Phys. B 320, 1 (1989); 
D. M. Asner et al. [CL EQ Collaboration], Phys. Rev. D 61, 012002 (2000) 
|arXiv:hep-ex/9902022|; E. M. Aitala et al. [E791 Collaboration], Phys. Rev. 
Lett. 86, 770 (2001) |arXiv:hep-ex/0007028| ; D. V. Bugg, Phys. Lett. B 572, 
1 (2003) [Erratum-ibid. B 595, 556 ( 2004)]; M. Ablikim et al. [BES Collabora- 
tion], Phys. Lett. B 598, 149 (2004) |arXiv:hep-ex/0406038| . 

[3] E. M. A itala et al. [E791 C ollaboration], Phys. Rev. Lett. 89, 121801 
(2002) larXiv:hep-ex/0204018| ; D. V. Bugg, Phys. Lett. B 572, 1 (2003) 
[Erratum-ibid. B 595 , 556 (2004)]; M. Ablikim et al. [BES Collaboration], 
|arXiv:hep-ex/0506055| 

[4] E. M. Aitala et al. [E791 Collaboration], arXiv:hep-ex705 07099 ; J. M. Link et al . 
[FOCUS Collaboration], Phys. Lett. B 621, 72 (2005) |arXiv:hep-ex/0503043| . 

[5] B. Aubert et al. [BABA R Collaboration], Phys. Rev. Lett. 90, 242001 (2003) 
|arXiv:hep-ex/0304021|; D. Besson et al. [ CLEQ Collaboration], Phys. Rev. D 
68, 032002 (2003) | arXiv:hep-ex/0305100| ; P. Krokovny et al. [B elle Collabo- 
ration], Phys. Rev. Lett. 91, 262002 (2003) | a rXiv:hep-ex/0308019| . 

[6] M. Nielsen, R. D. Matheus, F. S. Navarra, M. E. Bracco, and A. Lozea, 
|arXiv:hep-ph/0509131[ T. Mehen and R. P. Springer, Phys. Rev. D 72, 
034006 (2005) |arXiv:hep-ph/ 0503134) ; K. Terasaki and B. H J. McKel- 
lar, Prog. Theor. Phys. 114, 205 (200 5) |arXiv:hep-ph/050ll88 |. I. W. Lee, 
T. Lee, D. P. Min, and B. Y. Park, |arXiv:hep-ph/0412210[ C. Quigg, J. 
Phys. Conf. Ser. 9, 1 (2005) |arXiv:hep-ph/0411058|; "P~S. H wang and 
D. W. Kim, Phys. Lett. B 606, 116 (2005) larXiv:hep-ph/0410301|; Y. A. Si - 
monov and J. A. Tjon, Phys. Rev. D 70, 114013 (2004) [arXiv:hep-ph/040936l| ; 



7 



M. A. Nowak and J. Wasiluk, Acta Phys . Polon. B 35, 3021 (2004); 
A. Zhang, Phys. Rev. D 72, 017902 (2005) | arXiv:hep-ph/0408124|; Y. Q. Chen 
and X. Q. Li, Phys. Rev. Lett. 93, 232001 (2004) |arXiv:hep-ph/0407062| ; 
D. Becirevic, S. Fajfe r and S. Prelovsek, Phys. Lett. B 599, 55 (2004) 
|arXiv:hep-ph /0406296 ; P. Bicudo, Nucl. Phys. A 748, 537 (2005) 
|arXiv:he p-ph/0 401106|; E. E. Kolomeitse v and M. F. M. Lutz, Phys. Lett. 
B 582, 39 (2004) |arXiv:hep-ph/0307133"l; C. H. Chen and H. n. Li, Phys. 
Rev. D 69, 054002 (2004) |arXiv:hep-ph/0307075| ; Y. B. Dai, C. S. Huang , 
C. Liu and S. L. Zhu, Phys. Rev. D 68, 1140 11 (2003) |arXiv:hep-ph/0306274| ; 
G. S. Bali, Phys. Rev. D 68, 071501 (2003) |arXiv:hep-ph/0305209|; P. Colan - 
gelo and F. De Fazio, Phys. Lett. B 570 180 (2003) |arXiv:hep-ph/0305140| ; 
S. Godfrey, Phys. Lett. B 568, 254 (2003) |arXiv:hep-ph/0305122|; H. Y. Cheng 
and W. S. Hou, Phys. Lett. B 566, 193 (2003) |arXiv:hep-ph/0305038| ; 
T. Barnes, F. E. Close , and H. J. Lipkin, Phys. Rev. D 68, 054006 (2003) 
|arXiv:hep-p"h /0305025] . 

[7] E. van Beveren and G. Rupp, Phys. Rev. Lett. 91, 012003 (2003) 
|arXiv:hep-ph/03050"35] . 

[8] E. van Beveren and G. Rupp, Eur. Phys. J. C 22, 493 (2001) 
|arXiv:hep-ex/0106077| . 

[9] K. Abe et al [Belle Collaboration], Phys. Rev. D 69, 112002 (2004) 
|arXiv:hep-ex/0307021|; J. M. Link et a l. [FOCUS Collaboration], Phys. Lett. 
B 586, 11 (2004) |arXiv:hep-ex/0312060| . 

[10] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 94, 182002 (2005) 
|arXiv:hep-ex/0408126| . 

[11] E. van Bever en and G. Rupp, Mod. Phys. Lett. A 19, 1949 (2004) 
|arXiv:hep-ph /0406242] . 

[12] F. Kleefeld, AIP Conf. Proc. 717, 332 (2004) j|arXiv:hep-ph/0310320| . 

[13] E. van Beveren, G. Rupp, T. A. Rijken, and C. Dullemond, Phys. Rev. D 27, 
1527 (1983). 

[14] E. van Beveren an d G. Rupp, Phys. Rev. Lett. 93, 202001 (2004) 
|arXiv:hep-ph/04072"8T] . 

[15] http:/ /wwwasdoc.web.cern.ch/wwwasdoc/minuit/ 
minmain.html 

[16] The Mathematica Book, 5th edition, Wolfram Media, 2003, ISBN 1-57955-022- 
3. 

[17] A. V. Anisovich and A. V. Sarantsev, Phys. Lett. B 413, 137 (1997) 
|arXiv:hep-ph/970540T] . 

[18] E. van Beveren, C. Dullemond, and G. Rupp, Phys. Rev. D 21, 772 (1980) 
[Erratum-ibid. D 22, 787 (1980)]. 



8 



[19] T. Barnes, J. Phys. Conf. Ser. 9, 127 (2005) ar Xiv:hep-ph/0412057| . 

[20] Y. S. Kalashnikova, Phys. Rev. D 72, 034010 (2005) ||]arXiv:hep-ph/0506270| . 

[21] G. Rupp, F. Kleefeld, and E. van Beveren, AIP Conf. Proc. 756, 360 (2005) 
|arXiv:hep-ph/0412078| . 

[22] E. van Beveren, T. A. Rijken, K. Metzger, C. Dullemond, G. Rupp, and 
J. E. Ribeiro, Z. Phys. C 30, 615 (1986); E. van Beveren and G. Rupp, Eur. 
Phys. J. C 10, 469 (1999) |arXiv:hep-ph/9806246| . 

[23] R. Delbourgo and M. D. Scadron, Mod. Phys. Lett. A 10, 251 (1995) 
|arXiv:hep-ph/9910242] . 



9 



